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By Frederick C. Grant 


SLMMARY 


Linearized- theory calculations of the drag reduction achieved by 
applying the first three terms of a power series for twist to flat delta 
wings are presented. In addition, the reductions due to applying linear 
twist to a family of flat arrow wings are presented. The results cover 
the speed range of subsonic leading edges. 

The results show a 6-percent drag reduction due to twisting a flat 
delta wing with sonic leading edges and a steady decrease in the gains 
as sweepback increases. 

For the family of linearly twisted arrow wings investigated (that 
with sonic trailing edges), the maximum drag reduction is 2 percent in 
the medium sweepback range with a steady diminution in both directions. 

A beneficial effect of increasing aspect ratio obscures the twist effects 
in this case. The convergence to the optimum-power-series twist appears 
to be rapid. 


INTRODUCTION 


Numerous examples of drag reduction by warping sweptback supersonic 
wings have been calculated. (For example, see refs. 1 to 3-) In a recent 
paper (ref. 4), as many as 10 types of loading have been combined on a 
delta wing including power-series twist terms. In the present paper two 
cases of an arrow- wing plan form subjected to pure twist (spanwise slope 
variation) are considered. The twist is applied symmetrically from the 
root as a power of the distance from the root. Powers up to three are 
considered for a delta plan form. For an arrow plan form with sonic 
trailing edges the effects of linear twist are shown. 

The Lagrange method of reference 3 is used to compute both the 
optimum twist and the drag decrease over a flat wing for a given total 
lift. Leading-edge thrust is included in the drag calculations. 
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SYMBOLS 


b,f 

C D 

c D,ij 

( C D,ij)T 

C D,ij 

C L,i 

c L 

Cp 


E 

K 

k = \jl - k 2 
k = 0 q =' Pm 


proportionality constant 

parameters defining pressure distribution 

drag coefficient 

interference drag coefficient excluding leading-edge 
thrust 

interference drag coefficient of leading-edge thrust 


interference drag coefficient including leading-edge 
thrust 

lift coefficient of a wing with ith power twist 

lift coefficient of a combination of twisted wings 

part of lift coefficient due to ith-power component 
wing when Cl * 1 

lifting pressure coefficient (pressure on the lower 
surface minus pressure on the upper surface, divided 
by free-stream dynamic pressure) 

complete elliptic integral of second kind with modulus k 

complete elliptic integral of first kind with modulus k 


speed- sweepback parameter (see fig. l) 


M 


free-stream Mach number 


m 

S 

t = Py/x 
8 - py/x 
r\ = Py/x 


tangent of semiapex angle (see fig. l) 
wing area 


basic variables of generalized conical flow 
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V 

v n 

w 

x,y,z 

a 

(3 = \/m 2 - 1 
^C D 

1 - i* 

0 

Subscripts: 

f 

r 

t 

T 


free-stream velocity 

velocity of flow normal to leading edge, near leading 
edge 

velocity of flow in z-direction 

right-handed Cartesian coordinate system (see fig. l) 
local angle of attack 

flat-plate drag coefficient minus optimum drag coefficient 

root chord of arrow wing (see fig. 1) 
velocity potential 


flat 

power of twist variation 
dummy index or root 
tip 

plan form 


ANALYSIS 


Pressure Computation 

For a wing subjected to a symmetrical pure twist which varies as a 
power of the distance from the root (fig. 1), the downwash variation is 


w 

V 


-ai jyj 


i 


By introducing the variable T) = Py/x, equation (la) becomes 



(la) 


(lb) 
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The down-wash variation (eq. (lb)) is of the type treated in references 5 
and 6, and the general formula for the pressure distribution iB given in 
appendix A. If i is no greater than 3, the lifting pressure coeffi- 
cient is 


Cp,i 


/xY*' b i0 + b i2 e2 + b i4^ 

P \[ e o 2 - 92 


(6 = £y/x) (2) 


t 


The coefficients b-j_ r in equation (2) can be determined from, the fol- 
lowing matrix equation: 


CB = A 


(5) 


A general case of equation ( 3 ) is discussed in appendix A. For the case 
at hand, for linear twist. 



l( -1,2,0) 

l( -1,2,2) 


b 10 


4 

an 

P 

Cl- 

1 ( 0 , 2 , 0 ) 

1(0, 2, 2) 

B l = 

b 12 

A ! * 

0 

— 


for quadratic twist. 



- 


- 


1- 

i( -1,5,0) 

I(-1,5,2) 

b 2 = 

b 20 

Ag = 

- - a 2 
P 2 

I(l, 5,0) 

1 ( 1 , 5, 2 ) 

_ 


b 22 


0 





NACA TN ^4-104 


5 


and for cubic twist. 



— 


~ 


r i 


- 


l(-l,4,o) 

I( -1,4,2) 

l(-l,4,4) 


o 

IT\ 

n 


24 

a 3 

P 

C 3 = 

l(0,4,0) 

1(0, 4, 2) 

1(0, 4,4) 

b 5 = 

CVI 

fA 

.Q 

a 5 = 

0 


1(2,4, 0) 

1(2, 4, 2) 

1(2,4, 4) 


_ b 54_ 


0 


(4c) 


The l(a,b,e) functions of equations (h) are evaluated (appendix A) 
from the following relations: 


l(a+l,b+l,c) 


(c - b - l)l(a,b,c) 


" e 0 


I(a,b,c) 

d9 0 


(5a) 


l(a-b,0,c) = Ja-b+c ~ ^a-b+c+2 


(5b) 


Jn 



t n dt 

/i - 1 2 jt 2 - e 0 2 


(5c) 


If, for convenience, the following new variables are introduced: 


k = 0 O 



a±t = a i m± 

-(l+i-r) _ b ir 

4mai(ii) 


( 6 ) 
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equation (2) may be rewritten as follows: 


C 


Pi i 


4ma it (i.' )x* 


f i0 * f 12^ + f i4^ 
/l - ^ 


(T) 


By evaluating the elements of equations (4) and substituting in terms 
of the new variables (eqs. (6)), the coefficients f^^ are obtained 

as given in table I. 


Drag Computation _ 

Integration over the plan form yields the following equations for 
normalized lift and interference drag coefficients (4ma^(i!) = l) : 


<L,i 


1 

S 



2(1 - p ) 1+1 

i + 2 




(1 - m{r) i+2 |/l - ^ 


dijr 


( 8 ) 


mC D,ij 



1 (1 - u ) 1+J+1 

2 i + J + 2 




(9) 


It is understood _in equations (8) and (9) that only even values of r 
occur. 
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The interference drag given by equation (9) does not include leading- 
edge thrust which is calculated separately in appendix fe where it is 
shown that 


m 


(°Dj i j) T = 


4(1 - t-0 


i + j + 2 '■ 


f i0 + f i2 + f i4)( f j0 + f j2 + f j4) ^ 


The total interference drag may be written as the sum of equations (9) 
and (10) 

mCb, ±j = mC D>1 j + m(c D ^ij) T (ll) 


Calculated values of the drag-lift coefficients e.y 


°D,iJ 

C L,i C L,j 


derived from equations (6) to (ll) are given in table II. The quanti- 
ties and the analysis given in reference 3 allow the calculation 

of the optimum partition of a given lift among the twisted wings of this 
report for least drag of the combination. The optimum partitions of a 
unit lift coefficient among the twisted wings of this report, as well as 
the drag of these partitions, are calculated by the method of reference 3 
and are given in table III. 


RESULTS 


Delta Wings 

The total drag of the optimum twisted delta wings, together with 
the total drag of the component loadings taken alone, is shown in fig- 
ure 2(a); all are at unit lift coefficient. Unit lift coefficient on 
the twisted wings taken alone is attained by increasing the tip angle 
of attack a,^, the root angle of attack remaining zero. These values 
n~i 1 include leading-edge thrust. A lower bound for the drag calculated 
as in reference J is shown for comparison. It is pointed out in ref- 
erence 4 that this bound is a poor approximation to the mi nimum in the 
lower sweepback range. 

The small overall improvement in delta-wing drag level indicated 
by figure 2(a) is made by lowering the flat-plate wave drag at the 
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expense of increasing the vortex drag. The vortex drag is a minimum for 
the flat delta wing and is equal to the total drag at k = 0. The net 
result of this process in shown in figure 2(b) for a wing with optimum 
linear twist, as well as for the wing with optimum combined linear, quad- 
ratic, and cubic twist which is plotted in figure 2(a). Figure 2(b) 
shows that the greater part of the small gains obtained are contributed 
by the linear twist throughout the speed range. This indicates a rapid 
convergence to the optimum-power -series twist. At k = 1.0 the six 
additional cambered surfaces combined with the flat delta wing in ref- 
erence it gave about another 4-percent reduction in the total drag, or 
a total of nearly 10 percent. 


Arrow Wings 

Figure 3 is analogous to figure 2, except that in figure 3 only 
linear twist is combined with the flat plate. The family of arrow wings 
considered in figure 3 have sonic trailing edges. 

Figure 3(a) indicates that the possible gains due to twist are small. 
The increase in aspect ratio with increasing k has a beneficial effect 
in reducing the total drag which obscures the effect of the twist. The 
small twist effect is shown more clearly in figure 3(b) with a maximum 
drag reduction of 2 percent. 


Optimum Settings 

Figure 4 shows the settings required on linearly twisted wings for 
maximum drag reduction} CLg is the angle of attack required for unit 

lift coefficient on a flat wing, whereas a^. and cc^ are the root and 

tip angles of attack on the wing with optimum linear twist, also at unit 
lift coefficient. 


If a 60° delta wing at M = ]J2 


k-i] 




and a lift coefficient of 


0.2 are considered, figure 4(a) indicates that about 2° of washout is 
required. The drag reduction of 2.4 percent in this example compares with 
3.3-percent drag reduction that could have been obtained if quadratic and 
cubic twist had been included (fig- 2(b)). 


CONCLUDING REMARKS 


*• 


Calculations of the maximum drag reduction achieved by applying the 
first three terms of a power series for twist have shown that the drag of 
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* flat delta wings with sonic leading edges can be reduced approximately 

6 percent. The possible gains steadily diminish as sweepback is increased. 
Similar drag-reduction calculations for a family of linearly twisted arrow 

* wings with sonic trailing edges show a maximum drag reduction of about 

2 percent in the medium sweepback range with steadily decreasing gains in 
both directions. In this case, however, the beneficial effects of high 
aspect ratio obscure the effects of twist in the lower sweepback range. 

The convergence to the optimum-power-series twist appears to be rapid. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., June 13, 1957- 
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APPENDIX A 

GENERAL FORMULAS FOR PRESSURE DISTRIBUTIONS 


Generalized Conical Wings 

For dovnvash distributions of the rather general type 



where P^ is a polynomial of degree I, Heaslet and Lomax in refer- 
ences 5 and 6 have derived the lifting pressure coefficient 



(A 2) 


for the case of symmetrical sveepback. (See fig. 1.) The coeffi- 
cients b ir are determined in reference 6 by the following equation: 



(A5) 


Equating coefficients of T] on both sides of equation (A 3 ) yields 
i + 1 linear relations for b lr . Since solutions are additive, single 
terms of the polynomial P^ will be considered separately. 

Although equation (A3) is of concise and elegant form, it is not 
one suitable for calculations or convenient for discussion of solutions. 
In order to start a systematic reduction of the finite-part integrals of 
equation (A 3 ) , the following notation is introduced: 


l(a,b,c) 



(A4) 
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The derivative in the integrand of equation (Ah) is a homogeneous 
function of degree c - b - 1. For this derivative, Euler's relation 
for homogeneous functions is 



(A5) 

Multiplying equation (A5) by t a ^l - t 2 and integrating yields 

l(a+l,b+l,c) = (c - b - l)l(a,b,c) - 9 q l(a,b,c) (A 6) 

O0Q 


From equation (A6) it is plain that any l(a,b,c) can be derived 
by b applications of equation (A6) to 


l(a-b,0,c) - Ja.-t+c - J a -b+c+2 


(AT) 


where 


J n = 



(AB) 


In these equations Jn is, in general, an elliptic integral and as shown 
in reference 8 satisfies the following relation: 


(n + l)J n+2 - n(l + + (n - l)ls?J n _ 2 = 0 (A9) 

For n near zero, J n may be readily put in a standard form. 

Jn may be derived by using equation (A9) . For n 


near zero: 


Any other 
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J 


-1 


1 at 

E 2 


J 0 = -K 


J L = 


2 


J 2 = -E 





(A10) 


where the modulus of K and of E is equal to yl - k . Inspection of 
equation (A9) shows that Jn is a homogeneous linear function of E 
and K for even values of n. For odd values of n, J n is a rational 
function of IE. Since the derivatives of E and K 


dE 

dK 


~(K - E) 
k 2 


dK = k*TC - E 
dk Be 2 



(All) 


are also homogeneous linear functions of K and E, then l(a,h,c) is 
a homogeneous linear function of K and E, when a - h + c is even. 


If = aj_Tp, then a linear set of equations for h^ r may he 
written as follows: 


a i = 


J. 


i+1 
^ ' 




4(iJ) U ^ 
r=0 


> (A12a) 


i+1 

0= ^h lr l(j,i+l,r) (j=0,l,2,...,i-l) 
r=0 
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or in matrix form as 


l(-l,i+l,0) l(-l,i+l,l) 

l(0,i+l,0) l(0,i+l,l) 


l(-l,i+l,i+l) l-bio 
l(0,i+l,i+l) bn 




0 


l(i-l,i+l,0) l(i-l,i+l,l) 


l(i-l,i+l,i+l) 


*>ii 


h ±A+A 


which may also he written as 


CB = A 


0 


(AlZb) 


(A12c) 


Since the single term under consideration is an odd or even function 
as i is odd or even, the pressure must he an odd or even function 
with i. Thus, every other element of the B matrix must he zero. The 
first element is zero when i is odd; the second element is zero when i 

is even. Consider now the C matrix. Equations (A7) , (A9) , and (A10) 

show that, for odd values of p = a - b + c, l(a-h,0,c) is a polyno mi a l 

in — for p < 0 and a polynomial in 0 q for p > 0. The polynomial 
0 O 

degrees are -p for p < 0 and p + 1 for p > 0. When the recursion 
formula (A6) is applied i + 1 times to l(j-i,0,r), in every row except 
the first, the value of I is zero when p is odd. Thus, every other 

I is identically zero, except in the first row. When p is even, I is, 

in general, not zero. 

If the solutions of equations (A12) are now considered, every other 
equation, excluding the first, disappears. This means that for odd i 
and odd downwash or for even i and even downwash the number of variables 
and equations is the same and a unique solution exists. Reference 5 
points out that the form of the pressure solution (eq. (A2)) is unchanged 
if the natural sign of the downwash polynomial is required to change 
for T] < 0. With respect to equations (A12) , even solutions may he 
obtained from odd i, or odd solutions from even i, by reversing the 
natural choice of zero hj_ r in the B matrix. Again considering equa- 
tions (A12) , for odd i and even downwash or for even i and odd down- 
wash, it is found that the number of variables and equations is again 
the same and a unique solution again exists. 
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For a polynomial term of degree v < i, the linear equations are 


i+1 

0= ^ b lr l(j,i+l,r) (j=-l,0,...i-v-2,i-v,...i-l) 

r=0 

► 

i+1 

** ■ f (nW C ' 1)i ' V I *ir*i-v-l,i + l,r) 

r=0 J 


(A15a) 


or as before 


CB = A 


(Aljb) 


where C and B are unchanged, hut 


A = 


0 


jW 1 


o 


- v):v: 


(Aljc) 


For v < i, there is again a unique solution. for each v _ but it is not 
possible to demand odd solutions for even v_ or even solutions for odd. v. 
Such a demand sets to zero the right-hand side of equation (A13a) and 
there are no solutions . 


The underlying reason that the present^ analysis will not yield solu- 
tions in every case is that the vth and succeeding derivatives of the 
downwash terms with forced sign reversals have singularities at the origin 
which must be accounted for in the inversion of an integral equation. 

This was pointed out in reference 9 "by Kafka, who shows that a logarithmic 
singularity in the pressure at q = 0 is required to produce the required 
downwash. Therefore the solution to equation (A2) is incomplete when 
v < i. In order to include all cases, expressions analogous to equa- 
tions (A6) , (AT) , and (A8) may be derived by an analysis similar to that 
given, and the coefficients of the logarithmic terms may be determined. 
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Twisted Wings 

The downwash variation for a symmetrica ll y twisted wing is 



and this case is treated hy equations (A12) . For linear twist 

C1B1 = A X 

where 


I(-1,2,0) 

I(-1,2,2) 


— “ 
^lO 



Bi = 


1(0, 2,0) 

1(0, 2,2) 


13 12 


Use of equations (A6) , (A7) , (A9) , and (A10) yields 

I(-1,2,0) = - -Ar 
2k 5 

l(-l,2,2) = - = 

k I 

1(0, 2,0)= 

1 - k 2 


1(0,2, 2) 



In a similar manner, for quadratic twist 


(Alii) 


(A15a> 


4 

“ — a-i 

p 1 

0 

(A15b) 


(A16) 


C2B2 — A2 


(A17a) 
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and 


C 2 ~ 


l(-l,3,0) I (- 1 , 3 , 2 ) 

I(l,3,0) 1(1, 3, 2 ) 


b 2 = 


3 20 


3 22 


A2 = 


where 


I( -1,3,0) = + (2k 2 - 4k ^ )k 


-4 4 
Fk 


If-1,3,2) = (-6 + IQS 2 - +($? - jit ) 

k 2 ^ 

(4 - 2k 2 )e + ( -3 + k^K 


K 


1(1,3, 0 ) 


I (1,3/2) = 


And for cubic twist 


and 


_ (l + k^E - 2k 2 K 


C 3 B 3 ~ A 3 


| ®2 


0 

(A17b) 


(A18) 


(A19a) 



I (-1,4,0) 

I (-1,4, 2) 

I(-1,4,4) 


*30 


24 

"T 8 ? 

c 3 = 

1(0,4, 0) 

1(0,4, 2) 

1(0, 4, 4) 

B 5 = 

OJ 

KA 

P 

A 3 

0 


1(2,4, 0) 

1(2, 4,2) 

1(2, 4,4) 




0 


— 


- 

— 

- — 




(A19b) 
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where 


l(-l,4,0) = - 2* 

2)P 

l(-l,4,2) = - ^ \ 

k 5 


(A20) 




12 it 

k 
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APPENDIX B 

CALCULATION OF LEADING-EDGE THRUST 


In order to sustain a finite downwash at the subsonic leading edge 
of a wing, an infinite discontinuity in the pressure is required at the 
edge. Such an edge experiences a Buction force which depends on the 
nature of the pressure singularity (ref. 6). When the singularity in 
velocity normal to the leading edge has the form 


'n 


V 


Gn 


(Bl) 


the drag coefficient due to suction is given as follows: 


C D = 




(B2) 


In these equations, G n is a con stant, X y, is the perpend icular distance 
from the leading edge, and p n = \jl - where ^ = Nhy^l + m^. For 

the twisted wings of this paper (setting 4m(i.')ctat = 1) 



f i0 + f 12^ 2 + f i4^ 

ill - i 2 


(B3) 


Integration of equation (B3) from the leading edge in the manner of 
reference 6 yields 


where 



f iO a i+l + 


2 4 

a ^i2 c7 i-l a 


(B4) 
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Since the following relation holds: 


a r = 


- a^ r - 1 


a or. 


r-2 


any cr r can he reduced finally to the case r = 1 or r = 0- 
Performing this reduction in equation (B4) and differentiating with 
respect to gives 


Gh _ lim frr v n,i\ _ a 1+1 ( f i0 + f i2 + f i4) 

Pn "* 0\J n V 4 


where 


(B5) 


S = tan 




Substitution of equation (B5) into equation (B2) yields 

I 2 2 

(1 - p)mC D#i = - ^ ( f i0 + f i2 + f i4) W 

In the superposition of tw r o solutions , v n = v n ^ ^ + v n ^ j in equation (B5) 
and the interference drag is given by 

I Z2 

(1 - lx)ni(c D# ij) T = - J i 7 j +” 2 (^10 + f 12 + f iO ( f JO + f j2 + f j 4 ) 

(BT) 

Equation (B7) is easily extended for quartic, quintic, and higher degrees 
of twist by adding values of f inside the parentheses as required. 
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TABLE I 

VALUES OP POR TWISTED WINGS 


k 

f 00 

f 10 

f 12 

f 20 

f 22 

f 30 

f 32 

■^34 

0 

1 

0.63662 

0 

0.25000 

0 

g. 

0 

• 

0 

0 

0 

.2 

.95193 

.54453 

.04604 

.18650 

.04215 

.04524 

.01880 

.00016 

.4 

.86907 

.42435 

.10613 

.12271 

.08350 

.02463 

.03253 

.00103 

.6 

.78348 

.33026 

.15318 

.08282 

.10773 

.01480 

.03700 

.00248 

.8 

.70518 

.26189 

.18736 

.05822 

.12080 

.00924 

.03766 

.00423 

1.0 

.65662 

.21221 

.21221 

.04244 

.12732 

.00606 

.03638 

.00606 


/ 
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TABIE n 

INTERFERENCE DRAG-LIFT COEFFICIENTS FOE TWISIED WINGS 


[ m = l] 


(a) Delta wings, n = 0 


k 

e 00 

H 

O 

W 

e 02 

e 03 

e n 

e 12 

e 13 

€ 22 

6 23 

e 33 

0 

0.1592 

0.1592 

0.1592 

0.1592 

0.1790 

0.1910 

0.1989 

0.2122 

0.2274 

0.2487 

.2 

.1784 

.1834 

.1871 

.1900 

.2171 

.2417 

.2615 

.2848 

• 3215 

.3744 

.4 

.2204 

.2333 

.2422 

.2503 

.2926 

.3388 

.3792 

.4213 

.4965 

.6075 

.6 

.2790 

.3009 

.3167 

.3301 

.3932 

.4672 

.5286 

.5993 

.7147 

.8840 

.8 

.3559 

.3900 

.4170 

.4407 

.5259 

.6370 

.7291 

.8322 

1.0014 

1.2473 

1.0 

.5000 

.5683 

.6250 

.6744 

.7958 

.9834 

1.1399 

1.3000 

1.5746 

1.9662 


(b) Arrow wings, n = k 


k 

e 00 

e 01 

e ll 

0 

0.1592 

0.1592 

O.179O 

.2 

.1421 

.1483 

.1892 

.4 

.1353 

.1467 

.2104 

.6 

.1259 

.1383 

.2123 

.8 

.1022 

.1116 

.1744 

1.0 

0 

0 

0 
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TABLE HI 


OPTIMUM COEFFICIENTS OF UPP AND DRAG FOR TOISTED WINCS 
[a = C L = l] 

(a) Delta wing with. 0, 1 5 2 , and 3rd degree twist 


. C L,0 

C L,1 

cl , 2 


C DjO 

1.0000 

^0.0001 

“0.0001 

0 

O.O796 

1.0797 

.2592 

-.5958 

.2567 

.0888 

2.0980 

-3.3817 

3.7405 

-1.4569 

.1074 

2.0803 

-3.0406 

3.2043 

-1.2440 

.1348 

I.9886 

-2.7857 

3.0310 

-1.2339 

.1708 

1.8384 

-2.0685 

2.1257 

-.8956 

.2341 


(h ) Delta wing with 0 end 1st degree twist 


Cl 3 0 

Cljl 

Cd,0 

1.0000 

0 

0.0796 

1.1744 

-.1744 

.0888 

1.2774 

-.2774 

.1084 

1.3114 

-.3114 

.1361 

1.3347 

-.3347 

.1722 

1.4292 

-.4292 

.2353 


(c) Arrow wing with 0 and 1st degree twist 


C L,1 


C L,0 


1.0000 

1.178a 

1.21B7 

1.2012 

1.1759 

t’x 


C D,0 


0.0796 

.0705 

.0664 

.0617 

.0505 

0 


“Exact calculation would show these values to be zero. The values of Cp, 
are in general less accurately computed than the values of Cj^q. 

b Value of X may be assigned arbitrarily. 
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(a) Delta wings. \i = 0. 

Figure 4.- Optimum incidence angles on linearly twisted wings. 

m = C L = C L ^ f = 1. 







